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1.  Introduction 

We  consider  an  experiment  to  compare  v  treatments  in 
b  blocks  of  size  k  each.  The  model  of  observations  is 
assumed  to  be  linear,  additive  and  homoscedastic.  The 
problem  is  to  determine  an  optimal  design,  that  is,  an 
allocation  of  treatments  to  blocks  that  is  best  in  some 
sense. 

The  criterion  of  optimality  depends  on  the  objective  of 
the  experiment.  If  the  objective  is  to  compare  all  treatments 
with  one  another,  then  Kiefer,  in  a  series  of  well  known 
publications,  and  other  authors  have  determined  families  of 
optimal  designs.  One  celebrated  result  is  that  a  BIB  design 
is  optimal  in  any  reasonable  sense  (Kiefer  (1975)). 

In  many  experimental  situations,  all  comparisons  between 
treatments  are  not  equally  important.  For  example,  one 
treatment  may  be  a  standard  treatment  or  a  control,  enjoying 
a  special  status.  The  rest  of  the  treatments  (test  treat¬ 
ments)  are  to  be  compared  with  the  control  -  comparisons 
among  test  treatments  not  being  of  much  importance.  Some 
optimal  designs  for  these  experiments  are  now  available  in  the 
literature.  To  get  started  in  the  area,  the  reader  may  look 
at  Pearce  (i960),  Bechhoffer  and  Tamhane  (1981),  Majumdar 
and  Notz  (1983),  Giovagnoli  and  Wynn  (1983),  Constantine  (1983) 
Hedayat  and  Majumdar  (1984),  Jacroux  (1984),  Notz  (1984). 

This  is  a  list  of  some,  though  not  all,  important  papers. 


More  generally,  the  treatments  may  be  divided  into  two 
sets,  say  G  and  H  containing  u  and  w  =  v  -  u  treat¬ 
ments  respectively.  All  comparisons  of  two  treatments,  one 
from  G  and  one  from  H  are  important,  but  comparisons  within 
G  or  within  H  are  of  less  consequence.  In  case  one  treat¬ 
ment  is  a  control,  we  take  u  *  1.  Similarly,  experiments 
could  be  conducted  to  compare  several  treatments 
with  two  or  more  standard  treatments.  For  example,  a  factory 
may  be  using  two  types  of  equipments  for  a  job.  The  object 
of  an  experiment  could  be  to  determine  whether  some  or  all 
of  three  types  of  new  equipments  are  better  than  the  existing 
ones.  In  this  situation  of  several  controls,  both  u  and  w 
are  2  or  more.  Farmers  often  use  several  varieties  at  the 
same  time,  since  each  of  them  could  be  suseptible  to  a  different 
disease.  Thus,  in  agricultural  experiments,  the  object  could 
be  to  determine  which  "package"  of  varieties  perform  better; 
or  whether  some  variety  from  one  package  should  be  replaced 
by  a  variety  from  another. 

In  this  article  we  consider  the  problem  of  finding  optimal 
designs  for  these  and  similar  situations.  Cox  (1971,  p.  238) 
suggests  that  if  there  are  several  controls,  then  each  of  them 
may  be  used  once  in  each  block.  The  rest  of  the  design  would 
be  a  BIB  design  in  the  treatments  from  the  other  set.  We 
shall  formulate  the  problem  mathematically  in  the  tradition  of 
Kiefer.  Our  criterion  will  be  the  sum  of  the  variances  of 
the  treatment  differences,  one  from  each  set.  This  criterion 


is  statistically  meaningful  in  this  context.  We  shall  derive 
several  methods  of  obtaining  optimal  designs  -  applicable 
under  different  types  of  experimental  situations. 

In  Section  2  we  look  at  experimental  situations  with  blocks 
of  small  size;  that  is,  k  is  small  compared  to  u  and  w. 

We  develop  a  method  of  obtaining  optimal  designs  in  Theorem 
2.1  and  obtain  some  additional  refinements  for  the  case 
u  =  w,  in  Corollary  2.2.  When  u  =  w,  intuition  suggests 
that  the  best  design  might  have  half  of  each  block  filled  by 
treatments  from  G,  and  the  other  half  from  H  (some  adjustments 
being  needed  for  odd  k's).  While  we  can  essentially  prove 
this  for  k  2  4,  surprisingly  it  is  not  true  when  k  «  2. 

We  also  give  a  catalog  of  optimal  designs  for  some  values 
of  u  and  w  when  k  =  2  and  b  s  100. 

In  Section  3,  we  look  at  blocks  of  large  size.  The  method 
of  proving  optimality  is  entirely  different  here.  We  prove  a 
theorem  applicable  to  any  linear  model,  and  derive  optimal 
block  designs  in  a  corollary.  This  also  yields  previously 
unknown  results  for  the  case  of  a  single  control. 

The  general  theorem  of  Section  3  is  applied  to  models 
eliminating  two  sources  of  heterogeneity  in  Section  4.  This 
gives  us  optimal  row  -column  designs.  As  a  consequence  we 
could  generalize  some  results  of  Notz  (1984)  for  the  case  of 
a  single  control. 


We  consider  the  problem  of  comparing  v  treatments  in  b 
blocks  of  size  k  each.  The  treatments  can  be  divided  into 


two  sets  G  and  H  of  u  and  w  (=v-u)  treatments  respectively 
For  this  section  only,  we  assume  that  k  is  small;  indeed 
we  insist  that 

k  s  min  (u,  w)  (2.1) 

The  model  is  additive,homoseedastic  and  linear.  This  means 
that  if  Y^^  is  the  observation  on  treatment  i  (1  s  i  s  v) 
in  block  j  (1  *  j  *  b)  and  plot  4  (  1  s  K  *  k),  then 


Y 


ijt  “ 


+  0 


3  +  *«»' 


(2.2) 


where  are  assumed  to  be  uncorrelated  random  variables 

2 

with  mean  0  and  common  variance  a  .  The  unknown  constants 
and  represent  the  effect  of  treatment  i  and  the 

effect  of  block  j  respectively.  Let  £(u,  w,  b ,  k)  be 
the  set  of  all  possible  experimental  designs.  The  purpose  of 
the  experiment  is  to  compare  each  treatment  from  G  with 
each  treatment  from  H.  Comparisons  within  G  and  H  are  of 
secondary  importance  and  will  play  no  role  in  the  selection 
of  a  design.  We  want  to  choose  an  experimental  design  from 
S(u,  w,  b,  k)  which  minimizes 


£  Z 

g€G  h€H 


Var(T 


dg 


(2.3) 


as  d  varies  over  all  of  £(u  ,  w, b,  k),  Here 

denotes  the  BLUE  of  r  -  t.  under  design  d.  This  criterion 

g  h 

can  be  called  the  A-criterion  since  the  optimal  design 


minimizes  the  average  variance  of  the  contrasts  of  interest, 
tg  -  t^,  g  6  G,  h  €  H.  The  A-criterion  is  expected  to  be 
generally  acceptable  in  most  situations  since  it  is  statistically 
meaningful. 

If  d  €  £(u,  w,  b,  k),  then  we  denote  by  *  ndl^  *  k) 

the  number  of  times  treatment  i  occurs  in  block  j. 

Nd  =  (ndij^  is  the  incidence  matrix  of  d  and 
Cd  =  Ddr  -  (lA)NdN'  is  the  "c -matrix."  Here  Ddr  is  a 
diagonal  matrix  consisting  of  the  replications  of  the  treat¬ 
ments. 

Let  us  denote  by  1  and  I  the  vector  of  n  ones  and 

n  n 

the  n  *  n  identity  matrix  respectively,  and  let  »  1m1n* 

Write 

pi  - 1° . k . 0  I 

a  w  x  v  matrix  in  partitioned  form. 

p'  -  [Pi  ... 

then  (2.3)  can  be  written  as  a  function  of  Cd  in  the  form 

«(cd)  -  tr  P  C-P'  »  tr  P^-Pi, 

where  Cd  is  any  generalized  inverse  of  Cd.  The  following  is 
a  useful  subclass  of  vxv  matrices. 

C  ■  (C  :  C  is  nonnegative  definite,  rank  C  »  v  -  1, 

C  lv  *  0  and  for  some  scalars  p,  q,  r,  s  and  t,  C 


H 


can  be  written  as 


For  each  d  €  £(u,  w,  b,  k),  we  define 

b 

TdiJ  3  g|G  ndgj*  T<323  =  h6H  ndh«J*  Tdl=j=lTdir 

Td2  -  JjW  Sdl  -  Jj.  TdlJ  Sd2  “  J.\  TdSJ-  We  are 

now  in  a  position  to  state  and  prove  me  main  theorem  of  this 
section. 

Theorem  2.1.  Suppose  dQ  €  £(u,  w,  b,  k)  is  such  that 

(i)  C,  €  C,  (ii)  n  ^ ,  *  0  or  1  for  all  i,j  and 
a0  a0  J 

(iii)  f(Td  x  ,Sd  -l)  »  Min  f(Tdi,Sdl)  where 

aCT  a0  d€fl(u,  w,  b,  k) 

f(Tdrsdi)  58  ^(“di  “  sdi}  +  (w-i)24w(^-i)-k)Td2 
+  Sd2]  +  (u-l)2/[(u(k-l)-k)Tdl  +  Sdl],  then  dQ  is  A -optimal 

in  fl(u,  w,  b,  k) . 

Proof  Let  TT  be  the  set  of  all  permutations  n  of  the 
treatments  which  can  be  expressed  as  a  product  h  «  *r 

D  i  * 

where  rr  is  a  permutation  of  G  and  tv  a  permutation  of 
6  n 

H.  Let  Q(tr)  be  the  corresponding  permutation  matrix.  Then 
it  is  easy  to  see  that 


«(C,) 


Moreover,  0(0^)  is  convex  in  Cd.  These  two  properties  imply 
that 


*(C„) 


in  TT. 


u(u-l),  r  +  s 


where  E  is  over  all  rr 

u 

u 

A  °dgg/u'  «  - 

Z  £  c , 

g4g#  *1  dgg 

V 

,  s  *  EE 

cdhh'/w<w'1) 

h4h/=u+l 

h=u+l  °dhh/ 

a  v 

and  t=  2  S  c/,„Vl/uw.  Here,  we  write  G=  [1,2, 

g=l.h=u+l  agn 

H  =  (u+l,...,v}.  To  evaluate  )  we  must  find  the  eigen** 

* 

values  and  eigenvectors  of  Cd.  The  eigenvalues  can  be  obtained 
from  Constantine  (1981).  In  fact,  we  can  express  Cd  in  its 
spectral  form  as 


cd  a  pEl  +  rE2  " 


where  =»  (*u  u  q  ^uu 


-  (o 


*.-»• ’  UWT 


0> 

0  J* 

,  _  (*2  JUU  -  ufuw\ 

^  '  ^uwJwu  11 


The  Moore -Penrose  inverse  of  C.  is 

d 


Cd+  =  P  E1  +  r"1E2  -(vt)"1^. 

tf(Cd)  =  trace  PCd+P#  =*  w(u-l)/p  +  u(w-l)/r-  1/t,  with 

2 


lence, 
* 


A  A  (n«  ‘  Tdi/u>^(“-1), 


p  -  Tdl/u  - 


8 


b  v  2 

r  =  Tag/w  -  <ndhJ  "  Td2J//v')  A(w-l)  , 

*  -  -  Jx  TdlJ  Td2J/uWk  "  -<radl-Sdl)/mk- 

Let  us  fix  Tail*  •  •  •  »Tdlb  and  try  to  reduce  <*(Cd). 
Since,  by  condition  (2.1),  Tdl^/u  *  k/u  *  1, 

jfl  g^,  ^ndgj  ‘  ^1J  /U)2  *  Tdl  "  Sdl/u* 


with  equality  whenever  ndgj  *  0  or  1>  for  all  g  and  J. 
Similarly,  since  T^g^/w  *  1, 


b 

r 


*  (ndhj  '  Td2j/u,a  *  Td2  -  Sd2/"’ 


j*l  h^u+1 

with  equality  whenever  ndh^  =0  or  1,  for  all  h  and  j. 
These  relations  yield  a  lower  bound  to  $(Cd)  for  fixed 
values  of  TdH»  *  •  •  *^dlt*  A^'fcer  some  simplifications,  we  get 

<*(Cd)  *  0(C* )  *  u  w  f(Tdl,Sdl). 


Hence  the  theorem. 


Remark  1,  We  can  write  f(Tdi»Sdi)  as  a  func-tion  °f  Tdi 
and  Sdl  only  since  Tdl  +  Td2  =  bk  and  kTdl  -  Sdl  =  k  Td2 
To  carry  out  the  minimization  of  ^(^dl'^dl^  ^  could 

be  useful  to  introduce  integers  bda(asO, 1, . . . ,k) ,  where 


b,  is  the  number  of  blocks  for  which  T,. .  =  a.  In  this 


notation  b  *  s  b.  ,  T., 
o=0  dct  dl 


T.  ab  .  and  S,,  =  Z  cTbA. 
o=0  da  dl  o=0  da 


The  problem  reduces  to  minimizing  f(Tdi»Sdi)  38  f*(bd0,  •  •  •  »bdk) 

over  integers  b  (a*  0,l,...,k),  where  each  bdQ  =  0,1, ...,b 
k 

and  *  b  =  b#  If  k  is  small  then  this  problem  is  not 
o=0  d“ 

too  time  consuming.  For  example,  when  k  =  3,  u  =  3»  w  =  4 
and  b  a  30,  the  optimal  values  ba  of  bda  are  bQ  =  0, 
b^  as  18,  b2  at  12  and  b^  «  0.  For  a  clearer  representation, 
we  shall  denote  the  treatments  in  G  by  numbers  (1,2,...) 
and  the  treatments  in  H  by  alphabets  (A,B,...).  The 
following  design,  with  columns  as  blocks,  is  A-optimal  in 
5(3,4,30,3). 


112112112111111  222222  333333 
233233233AAABBC  AAABBC  AAABBC 
AAABBBCCCBCDCDD  BCDCDD  BCDCDD 


One  attractive  feature  of  this  example  is  that  only  two  optimal 
ba's  are  positive.  A  study  of  f(Tdl  Sdl),  with  some  examples, 
shows  that  this  may  not  be  true,  in  general.  Thus,  there  is 
no  guarantee  that  the  treatments  from  H,  as  a  collection,  are 
spread  as  homogeneously  as  possible  over  the  blocks,  which 
often  happens  when  u  =  1,  the  single  control  case  (Majumdar 
and  Notz  (1984)). 

When  k  =  2,  the  problem  is  to  minimize  l/bdl 
+  (u-l)/(bdl  +  (2  u/(u-l))bd2)  +  (w-l)/(bdl  +  (2  w/(w-l))bd0), 


10. 


over  integers  in  the  region  b  <}0  +  b  dl  +  bd2  =  b  *  the 
optimal  bda  is  denoted  by  ba,  then  the  necessary  and 
sufficient  conditions  for  obtaining  an  optimal  design  by  the 
method  of  Theorem  2.1  are  bQ  ■  0(mod  w(w-l)/2),  b^^  *  0  (mod  uw) 
and  b2  *  0  (mod  u(u-l)/2).  We  shall  denote  by  a-^ta,,  the 
set  of  all  ag-ples  which  can  be  chosen  from  a  set  of  a^ 
treatments.  Also  let  d^  denote  m^  copies  of  w  z  2  in 
the  treatments  from  H,  d^  denote  m^  copies  of  u  2  2 
from  G  and  d2  denote  m^  copies  of  the  uw  pairs  of 
treatments  formed  by  matching  a  treatment  from  G  with  a 
treatment  from  H.  Here  bQ  «  m^  w(w-l)/2,  =  m2  uw  and 

bg  *  m^  u(u-l)/2.  When  the  necessary  and  sufficient  conditions 
are  satisfied  the  A-optimal  design  is 

d  =  dQ  U  dx  U  d2. 

We  searched  and  found  A-optimal  designs  for  k  *  2  in  the 
range  v  <  i  s  100  for  treatments  (u,w)  =  (2,2),  (2,3), 

(2,4),  (2,5) ,  (2,6),  (3 ,3),  (3>4).  The  catalog  of  A-optimal 
designs  produced  by  Theorem  2.1  for  these  values  of  the 
parameters  are  given  in  Table  1.  This  table  gives  the  optimal 
values  bQ,  b^  and  bg  for  combinations  of  u,  w  and  b. 

For  example,  when  u  =  2,  w  =  3,  and  b  *  22,  Table  1  shows  that 
bg  a  3,  b^  =  18  and  b2  =  1.  Therefore  an  A-optiraal  design 
in  c(2,3»22,2)  is 


*  *L  *  a  •  .  •  #  •  ,  *  a  .*,*»*  •  "  *  *  »  *  »  *  »  “  »  ' 

»  *  «  *  *•  .  •  I  '  *  •  ,*•  ,  «  I  ’  ,  «  *  •  *  •  "  •  *  •  * 


1111222111222111222AAB 

2AAAAAABBBBBBCCCCCCBCC 

Actually,  our  search  was  more  extensive.  It  included  all  (u,w, 
b)  values  with  2  <  u  <  6,  u  <  w  <6  and  v  *  b  *  100. 

For  some  of  these  (u,w)  combinations  no  optimal  design 
could  be  obtained  by  Theorem  2.1  for  b  <  100. 

Given  u,w  and  b.  Theorem  2.1  could  give  an  optimal 
design.  In  case  it  fails  to  do  so,  it  still  gives  a  lower 
bound  to  the  A-value,  $(Cd),  for  a  d  €  £(u,w,b).  Usually 
the  experimenter  has  a  fairly  good  idea  of  what  an  efficient 
design  looks  like.  By  comparing  the  A-value  of  a  design  with 
the  lower  bound  given  by  Theorem  2.1  he  will  be  able  to 
make  a  good  assessment  of  the  efficiency  of  his  design.  At 
the  very  least,  this  is  expected  to  be  the  practical  usefulness 
of  Theorem  2.1. 

The  combinatiorial  restrictions  imposed  on  dQ  by  the 
condition  Cd  €  C  are  quite  severe.  Construction  of  an 
A-optimal  design  will,  in  general,  be  a  difficult  task. 

In  case  w  =  u,  one  can  obtain  a  simplification  of  Theorem 
2.1.  This  is  given  in  Corollary  2.2. 


Corollary  2.2.  Suppose  w  *  u  and  k  s  4.  Let  dQ  €  £(u,u,b,k) 
be  such  that 


(1)  n,  ,  ,  »  0  or  1.  i  m  1 . v.  .1  *  1. 
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(ii)  Cr 


pi  +  qJ 
r  u  uu 

N  T 


t  Juu  pIu  +  qJuu 


t  Juu  \ 

+  qJ 


for  some  p,  q  and  t. 


(Ill)  Td  ^  »  k/2  when  k  is  even 


[k/2J  or  [k/2]  +  1  when  k  is  odd 


for  each  j  ■  l,...,b,  where  [m]  denotes  the  largest  integer 
not  exceeding  m.  Then  dQ  is  A-optimal  in  £(u,u,b,k). 

Proof.  Let  d  e  d(u,u,b,k).  Start  by  defining  Cd,  the 
average  of  Cd  over  all  permutations  in  IT,  as  in  the  proof 
of  Theorem  2.1.  Then  define 

c*  «  (°  I[1)  c*(°  Iu> 

cdl  \la  0  )  cdUu  0  J 

and  Cd  »  (Cd  +  Cdl)/2.  Observe  that  d(Cdl)  =  $(Cd).  By 
the  convexity  of 


d(Cd)  *  d(Cd)  *  #(Cd*). 


Clearly,  Cd 


pi  +  qJ  t  J 
*  u  M  uu  uu  , 

\  t  J  pi  -f  qJ 

uu  v  u*H  uu 


v  V 

j,  where  p  +  q  =  S  c  /v, 


q  = 


t  = 


U  V 

(  £  £  c  +  £  £  c  ' 

v  fitfg'-l  dgg'  h4h/=u+ldhh 

(  £  £  c.  .  )/vr. 

xg=l  iv=u+l  dsh' 


')/2u(u-l) 


** 


tf(Cd  )  =  2u(u-l)/p  -  1/t. 

Some  computations  and  an  appeal  to  condition  (2.1)  yield. 


p  st  (bk(k(u-l)-u)  +  s  d)/ku(u-l) 


with  equality  whenever  n 


'did 


0  or  1  for  all  i,  j.  Here 

.  Recall  that  t  =  -  E  T..,  .  T.0 ,Au2=  -v,Au2 
b  >1  dlJ  d2J  «  d 

where  we  write  7d  =  £  Tdlj  Td2j  •  Since  7d  -  (bk-S  )/2» 

•1®1 


Sd  =  Sdl  +  Sd2 


** 


*(Cd  )  *  ku^(l/7d  +  4(u-l)V(buk(k-l)  -27  )}. 


It  is  easy  to  see  that.  Max  7d  is  achieved  when 

when  k  is  even 
[k/2]  or  [k/2]  +  1  when  k  is  odd. 


TdU  - 


for  each  J  =  l,...,b.  In  either  case.  Max  7d  <  "bk2/4.  Let 
us  study  the  function 


a(7d)  *  1/Td  +  4(u-l)2/(buk(k-l)  -  27d). 
a(7d)  is  convex  in  7d  and 

M'YdJ/frVd  *  al^d^^d(buk(k-l)  -27d)2) 


O  A 

where  a^Y^)  *  8(u-l)Yd  -  (  buk(k-l)  -  2Yd)  .  Using  some 

O 

calculus  one  can  show  that,  for  each  u  s  2,  a1(bk  /4)  <  0 

whenever  k  *  4.  Hence  a(Yd)  is  a  minimum  when  y<j  is 

a  maximum.  This  establishes  Corollary  2.2. 

Table  1  shows  that  Corollary  2.2  need  not  hold  when  k  *  2. 

For  k  ■  3»  the  question  is  open. 

In  some  situations  an  experimenter  may  be  interested  only 

in  comparing  the  sets  of  treatments  G  and  H  each  taken  as 

a  whole.  In  agricultural  experiments  the  object  could  be  to 

determine  which  of  the  groups,  G  or  H,  results  in  the  best 

overall  yield.  Here  the  problem  could  be  to  compare  the  averages 

Z  t  /u  and  S  th/w.  So,  we  have  to  choose  a  d  in 

g€G  &  heH  n 

.  .  A  A 

£(u,w,b,k)  which  minimizes  Varl  Z  Tri  / u  -  2  f,.v,/w). 

v  geG  g  heH  an 

Along  the  lines  of  Theorem  2.1  it  can  be  shown  that  if 

C,,  e  c  and  T,  , k/2  if  k  is  even  and  T ,  , .  »  [k/2] 
dQ  dQlj  dQlj 

or  [k/2]  +  1  if  k  is  odd  for  all  j,  then  dQ  is  optimal. 

In  fact,  no  restriction  on  k,  like  condition  (2.1),  is  necessary 
for  this  result. 

3.  Optimal  Designs  In  Large  Blocks 

In  this  section  k  will  usually  be  large,  so  that 
condition  (2.1)  is  not  assumed  any  more.  We  shall  start  by 
proving  a  general  result  for  Gauss  Markov  models,  and  then 
derive  A-optimal  designs  for  particular  models. 

Let  &  be  the  set  of  all  designs.  If  d  s  i  then  the 


observations  Y ,  will  have  the  Gauss  Markov  model 
a 

E(Yd)  -  Xld91  +  X2d32,  Var(Yd)  .  ah. 

o 

Here  the  vectors  8111(1  the  scalar  a  are  the  unknown 

parameters.  We  will  assume  that  X^Xld  is  nonsingular  for 
all  d  €  £«  Define 

Cd  s*XidXld  ’XldX2d^X2dX2d)  X2dXld 

If  *  is  a  real  valued  function  defined  on  nonnegative  definite 
matrices  then  dQ  will  be  said  to  be  f  -  optimal  for 
estimating  if 

*(QCd0Q')  *  CdQ#)»  for  al3L  d  €  L- 

Theorem  5.1.  Let  t  have  the  property  that  $(A-B)  *  *(A) 
whenever  A,  B  and  A-B  are  nonnegative  definite.  Suppose 
dQ  s  &  satisfies 

X2d0ICld0(XU0Xld0)'1<i'’0  (5 


*  ♦«J(X'1dXld)'V) 


for  all  d  €  C. 


Then  dQ  is  -  optimal  for  estimating 


Proof.  Suppose  A  is  a  positive  definite  matrix,  B  is  a 
nonnegative  definite  matrix  and  C  =  A  -  B  is  also  nonnegative 
definite.  Let  Q  be  a  matrix  with  ft(Q')  c  m(C),  where  ft* 
denotes  the  column  span  of  a  matrix.  Then  it  can  be  shown  that 

QC  V  a  QA*V  if  Bk'W  =  0. 


Hence  for  every  d  € 


♦(QC‘  Q'  )  *  *(Q(X'  Xld  ) -1Q# )  by  (3.1), 

a0  0  0 

<  *(Q(X#ldXld)"V)  by  (3.2), 

*  *(QCdQ/), 

since  there  exists  a  generalized  inverse  Cd  for  which 
Cd  “  ^XidXld^_1  is  nonne6ative  definite.  This  establishes 
the  theorem. 

We  shall  apply  this  theorem  to  the  block  design  model 
given  in  (2.2).  Here  9^  =  (*]/ •  •  •  >%)  and  9g  »  (Pr...,Pb). 
For  a  design  d  s  £(u,w,b,k),  X'ldXld  =  Ddr,  the  diagonal 
matrix  showing  the  replications  r<ji»  •  •  •  »rdv  of  the  treatments 
X2dX2d  *  kIb  3  XldX2d  a  Nd*  the  treatment  block  incidence 
matrix.  Let  the  treatments  be  divided  into  two  sets  G  and 
H,  as  in  Section  2,  and  suppose  Q  is  the  matrix  P  defined 
in  Section  2,  $  =  trace.  Equation  (3.1)  is  satisfied  whenever 


Nd 


t 


for  some  diagonal  matrix  &.  To  find  a  design  dQ 


which 


satisfies  (3.2)  one  has  to  minimize 


u  v  _  . 

E  S  (j1-  +  r1-) 

g»l  to*u+l  'rdg  rdh' 

over  all  designs  d  s  fi(u,w,b,k).  In  (3.4)  we  are  actually 

/  A  * 

minimizing  tr  V(P9ldJ*but  for  a  model  with  no  block  effects 
while  (3.3)  is  similar  to  a  treatment  -  block 
orthogonality  condition.  Corollary  3.2  gives  one  class  of 
optimal  designs. 

Corollary  3.2.  Suppose  k  ■  O(mod  (u  +  ./uw))  and 
k  ■  O(mod  (w  +  Vuw)).  Let  dQ  €  fl(u,w,b,k)  satisfy 


n<JQgJ  ■  +  "dQhJ 


k/(w  +  Jmi), 


for  g  =  l,...,u,  h  =  u  +  1,...,  v,  J  ■  ltw, . . . ,b .  Then 
dQ  is  A-optimal  in  fl(u,w#b,k). 

Proof.  Consists  of  verification  of  (3.3)  and  a  straightforward 
solution  of  (3.4)  using  Lagrange's  multipliers. 


Remark.  In  particular  when  u  *  1,  we  have  the  case  of  a 
single  control  in  G  and  w  ■  v  -  1  test  treatments  in  H. 
Corollary  3.2  says  that  as  long  as 


for  each  i  =  2,...,v  and  j  =  1, then  the  resulting 
design  is  A-optimal  for  control-test  treatment  comparisons.  H 
the  control  is  denoted  by  1,  and  the  test  treatments  by 


2,...,v.  To  our  knowledge,  the  result  is  not  available  in  the 
literature. 


4.  Optimal  Row  Column  Designs . 

Suppose  v  treatments  are  to  be  compared  in  an  rxc 
array.  The  response  in  row  j  and  column  t  is 


ijt  58  Pj  +  ^t  +  Ti  +  eijt! 


if  treatment  i  is  used  in  cell  (j,t);  i  =  l,...,v,  j  *  1,... 
t  as  l,...,c.  Here  p.  is  a  row  effect,  -n  is  a  column 

<3  v 

effect,  is  a  treatment  effect  and  is  a  random  error 

As  usual,  s's  are  assumed  to  be  uncorrelated,  with  zero 

2 

expectation  and  variance  a  .  The  treatments  are  divided  into 
two  sets  G  and  H  of  u  and  w  treatments  respectively. 
4(u,w,r,c)  denotes  the  set  of  all  designs. 

We  shall  apply  Theorem  5.1  to  this  setup.  Take 
9^  —  ( ... and  ^2  ”  ( ■ # « , n^,...,n^).  As  in 
Section  3»  d0  s  £(u,w,r,c)  is  A-optimal  if 


for  diagonal  A-^  and  Ag,  and  if  dQ  minimizes  the  expression 
(3.4).  For  a  design  d,  we  have  used  to  denote  its 

treatment-row  incidence  matrix  and  its  treatment -column 

Incidence  matrix.  In  the  following  corollary  (A)^  denotes 
the  (i,j)th  entry  of  a  matrix  A. 

Corollary  4.1.  Suppose  r  *  O(mod(u  +  ./uw)),  r  ■  O(mod(w  +  yuw)), 
c  ■  0(mod(u  +  ,7uw)  and  c  ■  O(mod(w +  7uw)) .  Let  dQ  e  £(u,w,r,c) 
satisfy 

(Nid0)gJ  -  c/(u+y™).  (»ldo)hj  =  e/(w+y™), 

(4. 

(N2dQ)gt  =  r/  ( u  +  Jwi)  t  (N2dQJht  "  r/(w  +  7uw), 

for  g  =  1,...,  a,  h  »  u  +  1,...,  v,  j  =  l,...,r,  t  =  l,...,c. 
Then  dQ  is  A-optimal  in  £(u,w,r,c). 

This  corollary  can  be  proved  along  the  lines  of  Corollary 
3.2.  Let  us  describe  one  class  of  designs  which  satisfy 
(4,3).  Let  m  be  an  integer  and  consider  a  latin  square  of 
order  s  =  m(w+yuw).  Change  any  m  symbols  of  the  latin  square 
to  A,  any  other  m  symbols  to  B,  and  so  on.  A,B,.... 


denote 


denote  the  treatments 


the  treatments  in  H,  while  1,2,... 
in  G.  Change  m(^w/u)  symbols  to  1,  m(./w/u)  symbols 
to  2,  and  so  on.  Call  the  resulting  square  L^.  Obtain 


L.,rt  by  starting  with  the  same  or  a  different  square,  and  doing 
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similar  operations.  Get  squares  L^,  i  =  1,2, 


V 


j  =  l,2,...,q2,  in  this  fashion.  Form  the  sq^  x  sq2  array 
L  =  Then  I*  is  an  A-optimal  design  in  £(u,w, 

sq1#  sq2).  There  are  plenty  of  pairs  (u,w)  for  which  this 
method  can  be  used  -  for  example  u  =  4,  w  =  16. 

This  method  of  construction  is  a  development  on  a  method 
used  by  Notz  (1984)  for  control-treatment  comparisons.  Indeed, 
Corollaiy  4.1  with  u  *  1  generalizes  Corollary  2.1  of  Notz. 


Our  proof  is  also  different.  Corollary  2.1  of  Notz  states  that 


if  one  starts  with  a  latin  square  of  order  p  +  p  and 


2  2 

changes  the  symbols  p  +  1, ...,p  +  p  to  0  (the  control) 


op  p 

then  the  resulting  design  is  A-optimal  in  £(l,p  ,p  +p,p  +p) . 


We  can  go  a  little  further  than  this  result  of  Notz  for 
control-treatment  comparisons.  If  L^(i  *  1,2,...,  q^j 


j  =  l,2,...,q2)  are  (p  +p)  x  (p  +  p)  squares  obtained  by 


Corollary  2.1  of  Notz,  then  the  array  ((L^))  is  A-optimal 
in  C(l,p2,(p2  +  p)q1,  (P2  +  P)q2).  For  example,  let  p  =  2, 


q^  »  2  and  q2  =  3.  Then  this  method  will  give  a  patchwork 
of  six  6x6  squares  which  is  A-optimal  in  £(1,4,12,18). 
There  are  other  methods  of  obtaining  an  A-optimal  design  too. 
For  example,  the  layout. 


Co 

S11 

Hail 

S21 

Co 

S22 

S31 

s32 

Co 

is  A-optimal  in  £(1,4,12,18) .  Here  Co  is  a  4x4 

square  consisting  only  of  the  control. (i  *  1,2,3  >  J  =  1,2,3) 

are  4x4  latin  squares  consisting  of  test  treatments  only. 

L1  and  Lg  are  6x6  squares  obtained  by  Corollary  2.1 
of  Notz  (1984). 

The  optimal  designs  of  Corollaries  3.2  and  4.1  have  a 

certain  model  robustness  property.  For  instance,  designs 

specified  by  (4.3)  remain  A-optimal,  by  Corollary  3.2, 

in  case  row  effects  are  absent  (p  *  0,  j  *  1,  2,...,  r), 

J 

and  even  when  column  effects  are  also  absent  (pj  **  0,  j  *  1,2,.. 
r\j.  —  0,  t  —  1,2, ...,c). 

The  technique  of  finding  optimal  designs,  given  in  Theorem 
3.1,  can  be  easily  extended  to  additive  models  eliminating 
several  sources  of  heterogeneity.  This  would  generalize 
Corollary  4.1. 

Finally  we  would  like  to  point  out  that  all  the  A-optimal 
designs  given  in  this  article  are  also  optimal  under  several 
other  criteria.  One  of  them  is  MV  -  optimality.  An  MV 
optimal  design  minimizes 

Max  Var(?.  -  r..  ) 


g€G,  h€H 


'dg  dh' 
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